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Nomenclature 
Letter Symbols 
H the pressure head placed into the fluid, ft 
J the advance ratio based on body diameter, nondimensional 
K an inlet energy loss coefficient, nondimensional 
P the static pressure, 1b/£t? 
¥ the free static pressure corresponding to the given 
submergence dept, 1b/ft2 
F the fluid vapor pressure at its bulk temperature, efee” 
Tr the rotor tip radius, ft 
tp ; the maximum body radius of the vehicle, ft 
s/C the space-to-chord ratio, nondimensional 
U the average peripheral blade velocity based on mass flow, ft/sec 
vi the average fluid meridional velocity component based on mass 
flow, ft/sec 
AV the average fluid velocity change between stations © and @ 
based on mass flow, ft/sec 
Vo the average fluid peripheral velocity component based on 
mass flow, ft/sec 
Vv the average absolute fluid velocity based on mass flow, ft/sec 
Ve the vehicle forward velocity, ft/sec 


Greek Letters 


a the absolute fluid angle, radians 

B the relative fluid angle, radians 

re) the fluid density, slug/ft> 

oO the cavitation index, nondimensional 


: 
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ee 
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Subscripts 
Ree ese ff 
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R 


Nomenclature (Continued) 


indicates 


indicates 


indicates 


indicates 


station location 
a mean value 
pertaining to the rotor 


pertaining to the stator 
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INTRODUCTION 
The following outlines a method for determining the blade space-to- 
chord ratio to provide a prescribed blade cavitation index. It also 
indicates the hydraulic efficiency that would be obtained as a function of 
the advance ratio and the ingested mass flow of the pumpjet. This procedure 
has been added to the preliminary design process described in fa", which 
in the past, has assumed the hydraulic efficiency at some constant value 
typical of axial flow pumps. As a first check, the Akron propulsor of [1] 
is reviewed to see if the procedure can be applied with reasonable results. 
In addition, a series of runs were made for a body having a drag coefficient 
of 0.103 and at a number of advance ratios for both a preswirl and conventional 


pumpjet. 


ESTIMATES OF EFFICIENCY AND BLADE CHARACTERISTICS FOR A CONVENTIONAL PUMPJET 
The efficiency of a stage of an axial flow compressor is presented in [2] 
and shall be used as a basis to predict the efficiency of the pumpjet based 
on section characteristics of the rotor and stator blading at the mean 
radius, rae The conventional pumpjet geometry is decribed in Figure 1. In 
the preliminary design phase a very limited knowledge of the flow field 
characteristics and blade section parameters are known and, therefore, it is 
necessary to make some engineering approximations. 
It is first necessary to obtain the space-to-chord ratio, or its 
inverse solidity, needed to satisfy the limits of cavitation that are 
desired. A conservative estimate is made that the static pressure forward of 


the rotor blade tip, Po> is equal to the free~stream static pressure, P.. 


* 
References Page 20 
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On this basis, the rotor blade pressure coefficient at cavitation inception, 


Cc. » becomes 
R 
P. = P p. -F (v 2 
Ys Vv 00 Vv © 
C Be ees = —— = Op lm, . (1) 
R Pu ai 4 2 
2 2 


The value of (Wo/V,,) is the relative velocity ratio approaching the 
rotor blade at the (0.9 rp/rp) radius and is defined by Equation (6). The 
impulse momentum relation can be used to approximate the average differnce 


in static pressure across the blade as outlined in [1]: 


Te c n 8 
ee, fe eee — Bae x (2) 


A review of a number of blade pressure diagrams has indicated the average 
pressure difference across a blade chord closely approximates the difference 
in pressure from blade inlet to the minimum pressure point on the blade. 


On this basis, 


be pv? |¥ Che Pal | Mss 
2 
The value of AV, /V,, for any ingested mass flow can be evaluated from: 
R 

™ (_ )? a 

V 2 V AV 
ee (tee t+ al tk a 2. = 
Vv oo ros) oo ie 
co 
2g 


(4) 
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where all the quantities in the brackets are tabulated on the ov=put of 


the preliminary pumpjet design program (PPDP) of [3]. The value (AV, /V,) 


R 
is confined to an upper limit of (U/V,). 
The quantity U/V,, can be expressed as: 
= r 
Be gs 
a (0.9 7? > > (5) 
0 B 


where rp/tp is computed in the (PPDP) for each mass flow. The value of 
(Wo/V) can be estimated as: 


2 2 
m og ‘ vn, | 
= is a (6) 


The preceding provides the necessary data to determine (S/C) from 
Equation (3) that satisfy cavitation requirements. An upper limit of 1.8 
is placed on this quantity. 

The estimate of blade drag is obtained as outlined in [2] and 
consists of increments due to profile and secondary flow drag. The annulus 
wall drag is not included since this has already been considered in the 
shroud surface skin friction drag with was added to the required propulsor 
thrust as indicated in [1]. 


The rotor blade drag is then: 


C =C. +€C 5 (7) 
DR Dp De 
and 
(s 
( = blade profile drag = 0.020 - 0.004 S| : 
D Cc 
iy R 
2 
Cy = secondary flow drag = 0.018 Cy ‘ 


s R 


= oe yn aid eel a in Pa ae nthe ng en 
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The value of the rotor lift coefficient Cy can be evaluated as: 
R 


8 
Cc =2 (5 : (8) 
R 


The total pressure loss coefficient through the rotor blade row can 


now be determined as outlined in [2] where: 


2 
AP. C cos Bo 
o> =C¢ = (9) 
Rg 2 DR s R cose 
ew, m 
2 


loss 


The values of B, and Bo in Equation (9) are estimated for the rotor in 


the following manner: 


= Vv 
tan 8, = = — i 
foe) ve 
2 
Vv 
ee: Von 
tan B, ys ee = ’ 
oo oo 1.15 V 
“2 


and 


ao 
tan BF = (tan B. + tan 83) F 


It is assumed in the above that the meridional velocity accelerates by a 
| factor of 1.15 in passing through the rotor. The preceding provides an 
estimate of total pressure loss through the rotor blade row and now it 


is necessary to perform a similar analysis for the stator. 
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The stator blade drag and solidity can be similarly evaluated, where 

in this case the stator is located behind the rotor. The relation for 
the average pressure difference over the stator blade is: (The subscript 


@) denotes a station at the stator leading edge.) 


ay gis foe Re ua 
pe as ae) ae Gn) 
PV, Sv 00 00 
po uy 
2 
where it shall be assumed that: 
¥ % 
ee ae 
pm bee ; 
co ec 
Vo Vo 
ee 1) } one 
< ’ 
Ve Ve 
and 
a pe — Ae % 
— V Vv 
oo 5 Pah: pod 
Vv V V i 
fee] foo) foe} 


It is necessary to now relate Equation (10) to the cavitation index of the 


stator blade system. This can be approximated by writing the energy relation 


between far upstream of the pumpjet and the stator blade inlet: 


m v- vy,” 
mm "5 9, — 
f, © RF fy te og te é (11) 
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The quantity (P) can be subtracted from each side of Equation (11) 


and the relation rearranged to give: 


Vy ek 1% ii 
Oo, =C —!| + jJ——j -|——|) + || - (12) 
Ss bg v. V M [v° Vv 2 
Cc 
2g 
The quantity C. is now expressed as, 
S 
V 
G¢ =2 §] Neo] | A 
bg (6 S v, Vo 
and can be substituted into Equation (12) to give: 
Dea: ne: 92 : 
Vv Vv f 
Gp lash eee ee it 
C} S Ve Ve; [r= “" 2 V, Vy 
= 4 
L = 
(13) 


The cavitation index substituted in Equation (13) is that value calcu- 
lated in the program for the rotor at the mass flow coefficient considered 
and assumes no acceleration or deceleration of the flow as it approaches 
the rotor. An upper limit for the space-to-chord ratio of 2.0 has been 


imposed on the output. The stator blade drag can be estimated where: 


¢. #i& #¢ c (14) 
D ye ily 
where 
C, = 0.020 - 0.004 a P 
| P Ss 
and 
Cc, = 0.018 .* : ‘ 
S S 
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For the stator assume: 
i AV 
S "a. 
) i) 4 
The total pressure loss coefficient for the stator is then: 
AP cos a 
Fe ee ace ae (16) 
° Vv é sl S cass 
PY, ‘m 
Z 
where: 
1.1V 
Re 
Ve 
tana, = ——_ ; 
+ 49 
nach 
an 
tan Oe = 0.0 $ 
and 
tana = = tan a 
m 2 4 i 
The hydraulic efficiency is then computed as: 
SR + bo bis) 
Wo 
Ny = 1 — oa ae a . (17) 
aaa Ven 
we (a | 
g 


| 
| 
j 


i 
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The preceding shall be applied to the pumpjet designed in [1] to determine 
efficiency, lift coefficient, etc., as determined by this technique and 
/ compared to that obtained in the Akron exercise at the selected design 
point. 
Appendix (I) lists the input data for the Akron pumpjet analysis and 
the new program. The output is also listed and consists of an output 
identical to that of [3] plus the output which indicates computed hydraulic 
efficiency, the rotor and stator space-to-chord, and their respective lift 
coefficients. The use of the roter lift coefficient in conjunction with 
Figure 2, which is an empirical correlation of lift coefficient and 
achievable blade pressure coefficient, will aid in selecting the optimum 
rotor tip detecer. The estimated hydraulic efficiency for the Akron pumpjet 
in [1] was 0.89, which closely agrees with that computed and listed as EFF-TOT 
in Appendix (I) at a r,/r, of 0.414. The lift coefficient of the rotor and 
stator is listed in Appendix (I) as 0.34 and 1.02 respectively. The values 


computed in [1] for the rotor and stator near the 0.9 Cp are 0.3 and 1.0. 


ESTIMATES OF EFFICIENCY AND BLADE CHARACTERISTICS FOR A PRESWIRL PUMPJET 

As depicted in Figure 3, this pumpjet arrangement would consist of a 
stator system upstream of a rotor. The stator places swirl in the flow 
counter to rotor rotation. In the case of the stator vanes, since the 
velocity over their surfaces is relatively low, the blade solidity is a 
function of blade loading rather than cavitation performance. On this basis, 
the upper limit of lift coefficient for the stators is assumed at 1.4. The 
output from the PPDP provides data to compute AV, /V., for each mass flow 


S 
from: 
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3 (al? ("a ; 
AV 2 AV 1 
: B vo a ye | 
Wy NV. =< __.. (18) 
\ 
S . Tr 
(2) (0.9) 4 aI 
B 


An upper limit has been imposed on %, /V., so as not to permit it to exceed 
3 
the blade speed at 0.9 r,/tp- 
The relation for the lift coefficient permits computing (S/C)., assuming 


the relative velocity ratio over the stator blade surface (V,/V,,) is unity. 


yu 
L Vv 
if} = Si - 22 (19) 
Ss Vo Vo 
eae 
Nie See 


At the larger values of ingested mass the values of (S/C). will be 
greater than 2.0. However, an upper limit of 2.0 has been imposed on 
this quantity. The preceding provides data to compute the profile and 
secondary drag coefficients for the stator. The profile drag coefficient j 
for the preswirl stator blade row is lower than that estimated for the 
stator of the conventional pumpjet because it is in an assumed uniform 
flow field that enters the blade row at zero incidence. On this basis, 


the drag coefficient is lower as indicated by [4]. 


where: 


\ 
C, = 0.012 - 0.004 a ; 
P S 
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The angles a, and on are approximated for inclusion in Equation (16) 


2 


to compute the total pressure loss coefficient through the stator as: 


The rotor solidity shall be estimated assuming that the relative velocity 


near the blade tip is: 


\2 


5 
aoe 
V 
fo) 


+ 


The value of Vy IY. is the same as computed for insertion in Equation (19). 
4 
The energy equation from upstream of the preswirl vanes to the inlet 
of the rotor, assuming the same meridional velocity at both stations can 


be expressed: 


(21) 


By subtracting vapor pressire (P.) from both sides of Equation (21) and 
using the relation given by Equation (3) the value of (S/C), can be written 


as: 
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It is assumed that zero acceleration or deceleration of the flow, from 
upstream to the rotor inlet, occurs when computing (o,) for insertion 
in Equation (22). On this basis, the energy equation can be written 
between these two stations and, using the relations for the minimum 


pressure coefficient, the cavitation index at rotor inlet is computed 


as: 


a: te: 
W, 2 "6, ‘my 
R oo) 0 i) 


An upper limit was placed on the above (S/C), value of 1.8. The lift 


coefficient for both the stator and rotor is estimated based on the 


f 
| 
| 
f 


relation: 

/ S “ey 
(3 = 2 |= —— (for stator) ; 
L C V 

S S oo 

AV, 
S 8 
Cc. = 2 §] < (for rotor) ; ] 
R R 4 


where the respective (S/C) of stator and rotor are inserted. 


The values of 8, and BF can be estimated for evaluation of Equation (9) 


by the following: 


2 “o, **t 
0.9 -— + — 0.9 —-— 
J TR ve J rR 1 
tan B, = = ; tan B. = Fz ; tan BF _ 2 (tan 8, = 8.) 
7‘ % 
4 4 
v 1.25 v 


ai%= 
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| 
where: | 
2 
2 4 = % 2 V | 
, aif Tr 9, be 
7 te 0.9 7 +—] + | rn 
os) TS ie a 


The preceding values are used to compute the drag coefficient of the 


preswirl rotor as: j 


DR Dp Ds 
where 
Ss 
Cy = 0.020 - 0.004 a ‘ 
P R 
and 
Cy = 0.018 c 2 
s L R 


The preceding quantities are then used to compute the total pressure loss 


coefficient indicated by Equation (9). The hydraulic efficiency of this 


preswirl stage is expressed as: 


= 
V 
2 
A 
R S W, 
t= 1-4 ————— : (24) 
oe 
2 2 
‘. 7a 
2g) |W, 


4 


A flow diagram of the preliminary design procedure is outlined in Appendix 


(11) for both the conventional and preswirl pumpjet configurations. 


ae 


i RAID Liane LA AS Sl 
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A sample problem has been run for the preswirl pumpjet configuration | 
and the conventional pumpjet. The input data and a listing of the program 
and output are shown in Appendix (III) for both types of pumpjets at advance 
ratios of 0.675, 0.844, and 1.125. The power requirements that would be 
obtained assuming a hydraulic efficiency of 88 percent versus that obtained 
from the described methods are plotted in Figure 4 for the conventional 


pumpjet and Figure 5 for the preswirl pumpjet. It is apparent from these 


ee a Fe eT 


plots that the optimum tip diameter and the efficiency varies with advance 
ratio for both types of pumpjets. 

The preswirl and conventional pumpjet are approximately equal in 
efficiency for a given advance ratio. The efficiency of either type of 
pumpjet iis to decrease at low values of advance ratio. This loss in 
efficiency is associated with higher blade surface velocities due to the 


higher shaft speeds at the lower advance ratios. 


SUMMARY 

The results of the preliminary design procedure described within 
provides data indicating the hydraulic efficiency and power requirements 
of a pumpjet operating on a given body at a specified advance ratio as a 
function of ingested mass flow. The space-to-chord ratio and lift 
coefficient of the blading near the tip (0.9 r,/t,) is computed and can 
be used in selecting blade chord, if blade number is known. The lift 
coefficient has been empirically correlated with achievable minimum blade 
pressure coefficients as indicated by Figure 2. This correlation in 
combination with the program output provides a guide in selecting the 


mass flow and lift coefficient that will provide the desired cavitation 


performance. 


I 
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The results of the described preliminary design procedure indicates 
that the maximum efficiency and the associated mass flow and tip radius 
of the pumpjet changes as a function of advance ratio. This is contrasted 


to past preliminary design procedures which assumed the hydraulic efficiency 


constant as a function of ingested mass flow. 
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Appendix I 


AKRON input data 


N = 12 
R/RB V/Voo 
0.677 0.540 
0.692 0.700 
0.708 0.790 
0.724 0.852 
0.740 0.897 
0.760 0.938 
0.780 0.965 
0.800 0.982 
0.830 0.998 
0.860 : 1.000 
0.900 1.000 
0.940 1.000 
CDBB = 0.052 
CDS) = 0.063 
vV1l0V2 = 1.15 
ETAH = 0.89 
PHI = 0.175 
LORB = 1.00 
RHL = 0.677 
RH2 = 0.250 
XXK = 0.13 
PRCNT = 0.0 
CB = 0.3 
PHI2 = 0.262 
NN = 1.0 
ADRAT = 1.0 


= CN . 4 
™ 
&- ; UEGSE*%O 1Stsero OS6721°9O L2520°O v2lve* v 100i °0 £39504%°0 Cele 7eL ouluce) CUN%6*D 
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Appendix IT 


Conventional Pumpjet Functional Diagram 


| SECTION ONE CALCULATIONS 


May 3, 1979 
WSG: JEF:mmj 


¥ 


| READ GEOMETRIC AND REFERENCE | 


VELOCITY DATA, DRAG COEFFICIENTS 
AND PSNED ERAAIL Ey HYDRAULIC EFFICIENCY 


Sa es ee een 


| SPLINE FIT THE REFERENCE VELOCITY PROFILE | 


+ 


CALCULATE THE THRUST COEFFICIENT, 
MASS FLOW RATE COEFFICIENT, AND 
VARIOUS VELOCITY RATIOS AS A 
FUNCTION OF RADIAL DISTANCE 
AT THE REFERENCE STATION 


¥ 


| READ A VALUE OF LOSS COEFFICIENT | 


Y 


CALCULATE PROPULSIVE COEFFICIENT 
AS A FUNCTION OF REFERENCE 
FLOW AREA AND ASSUMED 
HYDRAULIC EFFICIENCY | 


¥ 


| WRITE OUTPUT DATA | 


Proceed to Section Two Calculations 
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Conventional Pumpjet Functional Diagram 


SECTION TWO CALCULATIONS | 


¥ 


| READ PRCNT 


| 
| 
| CALCULATE NEW VALUES OF 
PROPULSIVE COEFFICIENT AND 
MASS FLOW RATE COEFFICIENT 
FOR THE DESIGN INLET AREA RATIO | 
> | 
CALCULATE THE ROTOR TIP RATIOS | 


¥ 
| READ A VALUE OF ADVANCE RATIO | 
¥ 


CALCULATE CAVITATION INDEX AS A 
| FUNCTION OF ADVANCE RATIO FOR 


EACH ROTOR TIP RATIO 


+ 


WRITE OUTPUT DATA 


+ 


Proceed to Section Three Calculations 
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Conventional Pumpjet Functional Diagram 
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SECTION THREE CALCULATIONS 


came 


| 
11 Values Determined as a Function of Reference Flow Area and Given Advance Ratio | 


Rotor Calculations 
+ 


CALCULATE NEW ROTOR TIP | 
RAIIO AND PERIPHERAL 
BLADE VELOCITY RATIO |. 


¥ 


CALCULATE PERIPHERAL 
ABSOLUTE VELOCITY RATIO 
MAXIMUM VALUE EQUAL | 


TO PERIPHERAL BLADE 
| VELOCITY RATIO | 


+ 


CALCULATE RELATIVE 
VELOCITY RATIO AT 
ROTOR INLET 


¥ 


| CALCULATE NEW 
| CAVITATION INDEX 


+ 


| CALCULATE ROTOR SPACE/ 
CHORD RATIO MAXIMUM 
VALUE OF 1.8 | 
+ 


| CALCULATE LIFT COEFFICIENT | 
| FOR THE ROTOR | 


+ 


| CALCULATE ROTOR eae 


a | ~ 


<_—. 


Return to SECTION TWO CALCULATIONS 
and read new value of advance ratio 


| + 
Stator Calculations 


Y 


CALCULATE HASS FLOW RELATED | 
VELOCITY RATIO, PERIPHERAL 
ABSOLUTE VELOCITY RATIO, 

AND THE RELATIVE VELOCITY 
| RATIO AT THE STATOR INLET | 


Y 


SPACE/CHORD RATIO 


| MAXIMUM VALUE OF 2.0 


: 
| 
1 | | 
| 
| 
| 
| 


CALCULATE STATOR | 
| 


CALCULATE LIFT 
COEFFICIENT FOR THE 
STATOR 


CALCULATE STATOR 


EFFICIENCY AND TOTAL 
| EFFICIENCY 


¥ 


CALCULATE NEW PROPULSIVE | 
COEFFICIENT | 


+ 


| WRITE OUTPUT DATA | 


ee a 
+ 


NO» ALL ADVANCE 
RATIOS BEEN 


| HAVE 
| ~~ USED? 
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Preswirl Pumpjet Functional Diagram 


| SECTION THREE CALCULATIONS 


All Values Determined as a Function of Reference Flow Area and Given Advance Ratio 


¥ 


Stator Calculations Rotor Calculations 
Y y 


SET PERIPHERAL ABSOLUTE 
VELOCITY RATIO AT ROTOR 
INLET EQUAL TO STATOR VALUE 


| CALCULATE NEW ROTOR TIP RATIO 
AND ROTOR PERIPHERAL VELOCITY 


| RATIO | 


¥ 


| CALCULATE PERIPHERAL ABSOLUTE | 


VELOCITY RATIO MAXIMUM VALUE 
EQUAL TO PERIPHERAL VELOCITY 


RATIO 


‘¢ 


Y 


CALCULATE ROTOR RELATIVE VELOCITY 
COMPONENT AND THE RELATIVE 
VELOCITY RATIO AT THE ROTOR INLET 
¥ 


| CALCULATE NEW CAVITATION INDEX 


¥ 


SET RELATIVE VELOCITY RATIO 
AT STATOR INLET EQUAL TO 1.0 
AND SET STATOR LIFT 
COEFFICIENT UPPER LIMIT AT 1.4 


| CALCULATE ROTOR SPACE/CHORD 
| RATIO MAXIMUM VALUE OF 1.8 


¥ 


| CALCULATE LIFT COEFFICIENT 
| FOR THE ROTOR 


¥ 


CALCULATE ROTOR EFFICIENCY 
AND TOTAL EFFICIENCY 


‘ 
CALCULATE NEW PROPULSIVE 
COEFFICIENT 
Return to SECTION TWO i ee | * 
EALuEE tone aed Feat SO RAIOS BEC EX | WRITE OUTPUT DATA | 
new value of advance USED?. ~ 
ee ee 


END  «——YES— 
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Appendix IIT 


With Drag Coefficient CDBB = 0.103 


N= 21 

R/RB V/Voo 
0.1188 0.430 
0.150 0.450 
0.200 0.485 
0.250 0.528 
0.300 C2575 
0.350 0.628 
0.400 0.680 
0.450 0.730 
0.500 0.775 
0.550 0.818 
0.600 0.860 

CDBB = 0.103 

cDS = 0.003 

VilOV2 = a0 

ETAH — @) . 88 

PHI = 0.1396 

LORB = 0.662 

RH1 = 0.1188 

RH2 = OelL 

XXK = 0.08 

PRCNT = 0.0 

CB = 0.3 

PHI2 = 0.174 

NN = 3.0 

ADRAT = 0.675, 

0.8 
Let 
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